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M.S. Bhagwat,! A. Holl,^ A. Krassnigg,^ CD. Roberts,^'^ and P.C. Tandy^ 

^ Center for Nuclear Research, Department of Physics, 
Kent State University, Kent, Ohio 44242 U.S.A. 
"^Physics Division, Argonne National Laboratory, Argonne, IL 60439-4843 U.S.A. 
' Fachbereich Physik, Umversitdt Rostock, D-18051 Rostock, Germany 

Features of the dressed-quark-gluon vertex and their role in the gap and Bethe-Salpeter equations 
are explored. It is argued that quenched lattice data indicate the existence of net attraction in the 
colour-octet projection of the quark-antiquark scattering kernel. This attraction affects the unifor- 
mity with which solutions of truncated equations converge pointwise to solutions of the complete gap 
and vertex equations. For current-quark masses less than the scale set by dynamical chiral symmetry 
breaking, the dependence of the dressed-quark-gluon vertex on the current-quark mass is weak. The 
study employs a vertex model whose diagrammatic content is explicitly enumerable. That enables 
the systematic construction of a vertex-consistent Bethe-Salpeter kernel and thereby an exploration 
of the consequences for the strong interaction spectrum of attraction in the colour-octet channel. 
With rising current-quark mass the rainbow-ladder truncation is shown to provide an increasingly 
accurate estimate of a bound state's mass. Moreover, the calculated splitting between vector and 
pseudoscalar meson masses vanishes as the current-quark mass increases, which argues for the mass 
of the pseudoscalar partner of the T(15) to be above 9.4 GeV. The absence of colour-antitriplet 
diquarks from the strong interaction spectrum is contingent upon the net amount of attraction in 
the octet projected quark-antiquark scattering kernel. There is a window within which diquarks 
appear. The amount of attraction suggested by lattice results is outside this domain. 
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I. INTRODUCTION 

It is difficult to overestimate the importance of dy- 
namical chiral symmetry breaking (DCSB) in determin- 
ing features of the strong interaction spectrum. For ex- 
ample, this phenomenon, which can even be expressed 
in models that omit confinement [1-3], is responsible 
for the remarkably small value of the ratio of tt- and p- 
meson masses and, inseparable from this, the generation 
of large constituent-like masses for dressed-quarks. It fol- 
lows both: that a sum of constituent-like dressed-quark 
masses sets a baseline for the masses of all light hadrons, 
except the pseudoscalar mesons; and that pseudoscalar 
"meson cloud" contributions are essential to an under- 
standing of hadron observables, such as the masses [4, 5] 
and form factors [6, 7] of octet and decuplet baryons. 

That stated, our primary concern herein is the dressed- 
quark-gluon vertex, T'^[q]p). We are generally interested 
in its form, how that arises, and the ways it affects and 
is affected by strong interaction phenomena. 

This vertex is a key element in the gap equation^ 

S{p)-^ = Z2{ii-p + 7n'''^) 

+ Z,J^ g^D^,{p~q)^j^S{q)ri{q,p), (1) 



[i] We employ a Euclidean metric, with: {7^, 7;^} = 25fj^; 7^ = 7fi; 
and a-b = Y.t^i 



and this brings an immediate link with DCSB. The gap 
equation has long been used as a tool for developing in- 
sight into the origin of DCSB, and searching for a con- 
nection between this phenomenon and confinement [8]. 
Moreover, the vertex is essential to a valid description of 
bound states and therefore to realising and understand- 
ing Goldstone's theorem, in particular, and current con- 
servation in general. These are additional qualities whose 
elucidation will form a large part of our discussion. 

The gap equation consists of other elements: D^,y{k), 
the renormalised dressed-gluon propagator; m*^™, the A- 
dependent current-quark bare mass that appears in the 
Lagrangian; and :— d^q/ {2t:Y , which represents 
a translationally-invariant regularisation of the integral, 
with A the regularisation mass-scale. The quark-gluon- 
vertex and quark wave function renormalisation con- 
stants, Zi(C^, A^) and Z2(C^, A^) respectively, depend on 
the renormalisation point, the regularisation mass-scale 
and the gauge parameter. 

The gap equation's solution is the dressed-quark prop- 
agator. It takes the form 

S{p)-^ = i-i-pA{p\e) + B{p\(') 

and provides direct access to the gauge invariant vacuum 
quark condensate 

.A 

-{qq)l= Yim Z^{Q\k^)N,ivj, S\qX) . (3) 
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where is the renormaHsation constant for the scalar 
part of the quark self-energy, through which the current- 
quark bare- mass is related to the running mass: 

Z2(C^A2)m'^-(A) = Z4(C^A2)m(C), (4) 

tr£) identifies a trace over Dirac indices alone and the 
superscript "0" indicates the quantity was calculated in 
the chiral limit, which is unambiguously defined in an 

asymptotically free theory: 

Z2(C^A')m'^'"(A) = 0, A»C- (5) 

The vacuum quark condensate is a primary order param- 
eter for DCSB [9]. 

It is a longstanding prediction of Dyson-Schwinger 
equation (DSE) studies [8, 10, 11] that the two-point 
Schwinger functions which characterise the propagation 
of QCD's elementary excitations are strongly dressed at 
infrared length-scales, namely, fc^ < 2GeV^; and it has 
become apparent that this feature is materially impor- 
tant in explaining a wide range of hadron properties [12]. 

Such dressing is also a feature of the dressed-quark- 
gluon vertex, a three-point function, and it is certain 
that the infrared structure of this vertex has a big im- 
pact on properties of the gap equation's solution, such 
as: multiplicative renormalisability [13-15]; gauge co- 
variance [14, 16, 17], and the existence and realisation 
of confinement and DCSB [1-3, 18-20]. For example, re- 
lated vertex Ansdtze, which agree in the ultraviolet, can 
yield solutions for the drcsscd-quark propagator via the 
gap equation with completely different analytic proper- 
ties and incompatible conclusions on DCSB [21, 22]. 

Dyson-Schwinger equation predictions for the be- 
haviour of dressed-gluon [23-25] and dressed-quark prop- 
agators [26] have been confirmed in recent numerical sim- 
ulations of lattice-regularised QCD [27, 28]. Indeed, de- 
tailed study provides an understanding of the circum- 
stances in which pointwise agreement is obtained [29, 30] . 
This level of sophistication does not prevail with the 
dressed-quark-gluon vertex, however. Acquiring that is a 
realisable contemporary goal, and it is to aspects of this 
task that we address ourselves herein. 

We have organised our presentation as follows. In 
Sec. II, after outlining some general properties of the 
dressed-quark gluon vertex, we recapitulate on a non- 
perturbative DSE truncation scheme [31, 32] that has 
already enabled some systematic study of the connection 
between the dressed-quark-gluon vertex and the expres- 
sion of symmetries in strong interaction observables. In 
doing this we are led to propose an extension of earlier 
work, one which facilitates an exploration of the impact 
that aspects of the three-gluon vertex have on hadron 
phenomena. To amplify the illustrative efScacy of our 
analysis we introduce a simple model to describe the 
propagation of dressed-gluons [33] that reduces the rele- 
vant DSEs to a set of coupled algebraic equations which, 
notwithstanding their simplicity, exhibit characteristics 
essential to the strong interaction. 



In Sec. Ill we capitalise on the simplicity of our model 
and chronicle a range of qualitative features of the 
dressed-quark-gluon vertex and dressed-quark propaga- 
tor that are common to our model and QCD. Of particu- 
lar interest are the effects of net attraction in the colour- 
octet quark-antiquark scattering kernel which we are able 
to identify. We follow that in Sec. IV with an analy- 
sis of the Bethe-Salpeter equations which can be con- 
structed, consistent with the fully dressed-quark-gluon 
vertex, so that the Ward-Takahashi identities associated 
with strong interaction observables are automatically sat- 
isfied. This property is crucial to understanding hadron 
properties and interactions [34-38]. In addition, we de- 
scribe the evolution of pseudoscalar and vector meson 
masses with growing current-quark mass. One outcome 
of that is a quantitative assessment of the accuracy for 
meson masses of the widely used rainbow-ladder trun- 
cation, which we determine by a comparison with the 
masses obtained with all terms in the vertex and kernel 
retained. Section IV also contains an analysis of colour- 
antitriplet diquark correlations. Here, in particular, the 
role played by attraction in the colour-octet projection 
of the quark-antiquark scattering matrix is noteworthy. 

We close our presentation with a summary in Sec. V. 



II. DRESSED QUARK-GLUON VERTEX 

A. General features 

This three-point Schwinger function can be calculated 
in perturbation theory but, since we are interested in the 
role it plays in connection with confinement, DCSB and 
bound state properties, that is inadequate for our pur- 
poses: these phenomena are essentially nonperturbative. 

Instead we begin by observing that the dressed vertex 
can be written 

ir«(p,g) = |A«r^(p,(z) =: rr^(p,g); (6) 

viz., the colour structure factoriscs, and that twelve 
Lorentz invariant functions are required to completely 
specify the remaining Dirac- matrix- valued function; i.e., 

r^(P,a) = +7 • b + 'Z) (P + 9Vf2(p,g) 

-iip + q)f.h{p,q) + [■■■], (7) 

where the ellipsis denotes contributions from additional 
Dirac matrix structures that play no part herein. Since 
QCD is rcnormalisable, the bare amplitude associated 
with 7^ is the only function in Eq. (7) that exhibits an 
ultraviolet divergence at one-loop in perturbation theory. 

The requirement that QCD's action be invariant under 
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local colour gauge transformations entails'"! [39] 

J^B(e) {[1 - B{p, q)] S{p)-' - S{q)-' [1 - B{p, q)]} , 

(8) 

wherein T^{k!^), k = p — q, is the dressing function that 
appears in the renormalised covariant-gauge ghost prop- 
agator: 



fc2 



(9) 



and B{p,q) is the renormalised ghost-quark scattering 
kernel. At one-loop order on the domain in which per- 
turbation theory is a valid tool 



^«(fc^c') = 



7g//3l 



(10) 



with the anomalous dimensions 7g = — (3/8)C2(G) in 
Landau gauge, which we use throughout, and f3i = 
-(11/6) C2(G) + (1/3) A^/, where C2(G) = A^c and Nf 
is the number of active quark flavours. This result may 
be summarised as J^^{k'^) « 1 on the perturbative do- 
main, up to Inp^/AqQ^-corrcctions. In a similar sense, 

B{p,q) fa in Landau gauge on this domain. ['"1 

Equation (8) is a Slavnov- Taylor identity, one of a 
countable infinity of such relations in QCD, and it is 
plainly an extension of the Ward-Takahashi identity 
that applies to the fermion-photon vertex. The Ward- 
Takahashi identity entails that in the vertex describing 
the coupling of a photon to a dressed-quark: 



tup,p) = A{p\e). 

_d 
dp- 
_d 

dp' 



(11) 

(12) 
(13) 



Identifying and understanding this nontrivial structure 
of the dressed-quark-photon vertex has been crucial in 
describing the electromagnetic properties of mesons [40- 
43]. 

The similarity between the Slavnov- Taylor and Ward- 
Takahashi identities has immediate, important conse- 
quences. For example, if the result 



0<jrs(A;2)[i_B(p,g)]<oo 



(14) 



[ii] NB. Equation (8) is equally valid when expressed consistently in 
terms of bare Schwinger functions. 

[ill] An even closer analogy is a kindred result for Z(p^) in Eq. (2); 
viz., in Landau gauge [1 — Z(jp,(^^)] s at one loop in perturba- 
tion theory and hence, on the perturbative domain, corrections 
to this result arc modulated by Inlnp^/Aq^j-,. This very slow 
evolution is exhibited, e.g., in the numerical results of Ref [19]. 



also prevails on the nonperturbative domain then, be- 
cause of the known behaviour of the dressed-quark propa- 
gator, Eq. (14) is sufficient grounds for Eq. (8) to forecast 
that in the renormalised dressed-quark-gluon vertex 



1 < Ti{p,p) < oo 



(15) 



at infrared spacelike momenta. This result, an echo of 

Eq. (11), signals that the complete kernel in the DSE 
satisfied by (p, q) is attractive on the nonperturbative 
domain. 



B. Vertex in the gap equation 

The ability to use the gap equation to make robust 
statements about DCSB rests on the existence of a sys- 
tematic nonperturbative and chiral symmetry preserving 
truncation scheme. One such scheme was introduced in 
Ref. [31]. It may be described as a dressed- loop expansion 
of the dressed-quark-gluon vertex wherever it appears in 
the half amputated dressed-quark-antiquark (or -quark- 
quark) scattering matrix: S^K, a renormalisation-group 
invariant, where K is the dressed-quark-antiquark (or - 
quark-quark) scattering kernel. Thereafter, all n-point 
functions involved in connecting two particular quark- 
gluon vertices are fully dressed. 

The effect of this truncation in the gap equation, 
Eq. (1), is realised through the following representation 
of the dressed-quark-gluon vertex [32] 

r^{k,p) = 7^ + C^{k,p) + C+{k,p) + [...], (16) 



with 



xjpS{£ + k-ph^Sieh„, 



(17) 



^t{k,p) = ^J^ g^D,,^{£)DMe + k-p) 
X 7,, S(p - I) -fc, k-p), (18) 

wherein F^^ is the dressed-three-gluon vertex. It is ap- 
parent that the lowest order contribution to each term 
written explicitly is 0{g^). The ellipsis in Eq. (16) rep- 
resents terms whose leading contribution is 0{g^); e.g., 
crossed-box and two-rung dressed-gluon ladder diagrams, 
and also terms of higher leading-order. 

The C2 term in Eq. (16) only differs from a kindred 
term in QED by the colour factor. However, that factor 
is significant because it flips the sign of the interaction in 
this channel with respect to QED; i.e., since 



l<^l''l- = {lc2{G)-C2{R)\l'> = ^l' 



cf. l^lcl" = -C2iR)lc = - 



Nl-^ 
2Nc 



(19) 
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then single gluon exchange between a quark and anti- 
quark is repulsive in the colour-octet channel. Attrac- 
tion in the octet channel is provided by the £j term in 
Eq. (16), which involves the three-gluon vertex. These 
observations emphasise that Eq. (8) cannot be satisfied if 
the contribution from the three gluon vertex is neglected 
because the Slavnov- Taylor identity signals unambigu- 
ously that on the perturbative domain there is net at- 
traction in the octet channel. 

It is apparent, too, that the term involving the three 
gluon vertex is numerically amplified by a factor of N'^ 
cf. the C2 (Abelian-like) vertex correction. Hence, if the 
integrals are of similar magnitude then the TV^-enhanced 
three-gluon term must dominate in the octet channel. 
This expectation is borne out by the one-loop perturba- 
tive calculation of the two integrals exhibited in Eqs. (17), 
(18) and, moreover, the sum of both terms is precisely 
that combination necessary to satisfy the Slavnov- Taylor 
identity at this order [44]. 

C. Vertex model 

In illustrating features of the nonperturbative and 
symmetry preserving DSE truncation scheme introduced 
in Ref. [31] in connection, for example, with DCSB, con- 
finement and bound state structure, Ref. [32] employed 
a dressed-quark-gluon vertex obtained by resumming a 
subclass of diagrams based on C2 alone; namely, the ver- 
tex obtained as a solution of 

r;(fc+,fc_) = zr'j^ + ^j\'D,4p-e) 

xjpS{e+)T-{£+,e.)S{£-h.. (20) 

It was acknowledged that this subclass of diagrams is 
l/A^c-suppressed but, in the absence of nonperturbative 
information about C2 in general, and the dressed-three- 
gluon vertex in particular, this limitation was accepted. 

Herein we explore a model that qualitatively ame- 
liorates this defect while preserving characteristics that 
make calculations tractable and results transparent; viz., 
in Eq. (16) we write 

C2+Ct'^^2, (21) 

where 

.A 

£g(fc,p) :=-CC2(ii) / g^Dp,{p-e) 
Je 

and work with the vertex obtained as the solution of 

.A 

rC(fc+,fc_) = zrS^- cc2(i?) / g^Dp,{p-e) 

X 7,5(£+)r^(f+,£_)5(£_)7a. (23) 

To explain this model wc remark that the parameter 
C is a global coupling strength modifier. (NB. The value 



C = —(1/8) reproduces the vertex rcsummed in Ref. [32].) 
It is introduced so that our Ansatz may mimic the effects 
of attraction in the colour-octet channel without speci- 
fying a detailed form for the three-gluon vertex. This 
expedient will give a faithful model so long as the inte- 
grals over the momentum dependence of C2 and £ J that 
appear in our calculations are not too dissimilar. This is 
plausible because: they are both one-loop integrals pro- 
jected onto the same Dirac and Lorcntz structure and 
hence are pointwise similar at this order in perturbation 
theory; and their direct sum must conspire to give the 
simple momentum dependence in Eq. (8). Moreover, as 
we shall demonstrate, the model has material illustra- 
tive capacity and that alone is sufficient justification for 
proceeding. 

As we have already noted, the value C — —(1/8) cor- 
responds to completely neglecting the three-gluon vertex 
term. There is also another particular reference case; 
namely, C = 1. In this case a dressed-quark propagator 
obtained as the solution of the rainbow truncation of the 
gap equation: 

SR{p)-^^Z2{i-/-p + m}"^) 

+ ^ g^Di,^{p - q) y7^ SR{q) y7, , (24) 

when inserted in Eq. (23), yields a vertex r'^^{k,p) that 
satisfies 

{k - p)^ trl''{k,p) = SRik)-' - Snip)-' ; (25) 

viz., a Ward- Takaliashi identity. This is not materially 
useful, however, because herein we will seek and work 
with a dressed-quark propagator that is not a solution 
of Eq. (24) but rather a solution of Eq. (1) with a fully 
dressed vertex, and that vertex, determined self consis- 
tently, will not in general satisfy Eq. (25). 

D. Interaction model 

A simplification, important to our further analysis, 
is a confining model of the dressed- gluon interaction in 
Eq.(23). We use [33] 

V^^k) := g^D^^k) = - ^) {27rfg^ S\k) . 

(26) 

The constant G sets the model's mass-scale and hence- 
forth we mainly set ^ = 1 so that all mass-dimensioned 
quantities are measured in units ofQ. Furthermore, since 
the model is ultraviolet-finite, we will usually remove the 
regularisation mass-scale to infinity and set the renor- 
malisation constants equal to one. 

In these things we follow Ref. [32] . In addition, we 
reiterate that the model defined by Eq. (26) is a pre- 
cursor to an efficacious class of models that employ a 
renormalisation-group- improved effective interaction and 
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whose contemporary application is reviewed in Refs. [10- 
12]. It has many features in common with that class 
and, in addition, its distinctive momentum-dependence 
works to advantage in reducing integral equations to al- 
gebraic equations that preserve the character of the orig- 
inal. There is naturally a drawback: the simple mo- 
mentum dependence also leads to some model-dependent 
artefacts, but they are easily identified and hence not 
cause for concern. 



E. Algebraic vertex and gap equations 

If Eq. (26) is used in Eq. (23) then that part of the 
vertex which acts in the gap equation has no dependence 
on the total momentum of the quark- antiquark pair; i.e., 
only r^(p) := r^(p,p) contributes. In this case just the 
terms written explicitly in Eq. (7) are supported in our 
model for the dressed vertex, which can be expressed 



This is not a severe handicap because these Dirac struc- 
tures arc precisely those which dominate in Eq. (8) if 
B{p,q) is not materially enhanced nonperturbatively. 
The vertex equation is 

r^(p) = 7,. - C 7, Sip) r^(p) Sip) 7, , (28) 

where we have used the fact that C2{R) = 4/3 when 

Our analysis now mirrors that of Ref. [32]. The solu- 
tion of Eq. (28) is: 



(29) 



i=0 



= m [7m"m(p^) + 7-PPi^a2,iiP^)- iPi^Oi'i^iiP^)] , 

(30) 

where ^(p) satisfies a recursion relation: 

= -Cl^ Sip) Tl,{p) S{p) J, , (31) 
with r^^^^o = 7(ii) bare vertex, so that 

aio = 1 , a^.o = = a|o • (32) 

In concert with Eq. (30), Eq. (31) yields an algebraic 
matrix equation (s = p^) 

af+i(s):= I al,^,is) \ = [s; A, B) cx'i is) , (33) 
(4,i+ii.s) 



where (A(s) = sA^is) + B'^is) ) 
0^is;A,B) = 



2C 
■ A2 



/-AO ° \ 

2^2 5^2-^2 2AB 

\ AAB AsAB 2(^2 _ sA^) J 



(34) 



It follows from Eqs. (30) and (33) that 



1 



l-O^ 



(35) 



and hence, using Eq. (34), 
A 



af = 



Q!2 



A-2C' 



4CA2 



(A-4C) 



A2 + 2C(S2 - sA2) - 8C2 (A - 2C) 



SCAB 



as 



A2 + 2C(S2 - sA^) - 8C2 ■ 



(36) 
, (37) 



(38) 



With these equations one has a closed form for the 

dressed-quark-giuon vertex. 

It is evident that the momentum dependence of the 
vertex is completely determined by that of the dressed- 
quark propagator whose behaviour, however, the vertex 
itself influences because it appears in the gap equation: 

S(p)-^=i^-p + m + j^,Sip)Tlip). (39) 

Subject to Eq. (27), the gap equation expresses two cou- 
pled algebraic equations: 



Ais) = 1 + 



sA"^ + S2 



[Ai2a\-sa^)-Ba^], 



(40) 



^(*) = ""+ 5^2+ g2 (4^1 + ) - SA 4] . 



(41) 



The dressed-quark propagator and -quark-gluon vertex 
follow from the solution of these equations, which is gen- 
erally obtained numerically. 

In the chiral limit, which here is simply implemented 
by setting to = 0, a realisation of chiral symmetry in the 
Wigner-Weyl mode is always possible: it corresponds to 
the B = solution of the gap equation. However, since 
the phenomena of QCD are built on a Nambu-Goldstonc 
realisation of chiral symmetry, we do not consider the 
Wigner-Weyl mode any further. Its characterisation can 
be achieved in a straightforward manner by adapting the 
analysis of Ref. [32] to our improved vertex model. 
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III. PROPAGATOR AND VERTEX SOLUTIONS 
A. Algebraic results 

At this point some observations are useful in order to 
establish a context for our subsequent results. To begin 

wc explore the ultraviolet behaviour of the model. It is 
ultraviolet finite and hence at large spacelike momenta, 
s » 1 (in units of Q'^) 



A{s) 
Bis) 



1 + 



m + 



Si 

S 

bi 



(42) 
(43) 



with m the model's finite current-quark mass. The model 

is useful because these results persist in asymptotically 
free theories, up to Inp-^/Aq^p-corrections. With this 
behaviour it follows from Eqs. (36) - (43) that 



a?(5) 
a^(s) 



1 



2C 

s 



4C 



8C 



(44) 
(45) 
(46) 



and these results in turn mean that in the ultraviolet 
the behaviour of the massive dressed-quark propagator 
is determined by the ai term in the vertex, so that 



a,i = 2, bi = 4m . 



(47) 



The expansion of ai{s) around 1/s = reported in 
Eq. (44) is the same as that which arises in QCD, apart 
from the usual ln(s/AQ(~,j-,)-corrections. However, the 
leading terms in 02,3 are different: on the perturbative 
domain in QCD these functions both begin with a term 
of order (1 / s)[\n{s / Aq^j-,)]'^ , with d some combination of 
anomalous dimensions. The reason for the mismatch is 
readily understood. At one-loop order in QCD 



ln(//A|cD) 



7r//3i 



(48) 



with m the renormalisation point invariant current-quark 
mass and 7™ = (3/2)C2(i?). (This makes explicit the 
logarithmic correction to the leading term in Eq. (43).) 
For the purpose of this explanation then, on the per- 
turbative domain, with !F^{k^) « 1 and B{p,q) w 0, the 
Slavnov- Taylor identity, Eq. (8), is approximately equiva- 
lent to the Ward-Takahashi identity. Hence, via Eq. (13), 



"3(5) ~ — 

s fii 



-ln(pVA^CD) 



7r//3i 



(49) 



It is thus evident that in QCD, even though they are not 
themselves divergent, the leading order terms in both 

Q^2,3 are induced by the momentum-dependent renor- 
malisation of elements contributing to their evaluation. 



Such terms are naturally missing in our ultraviolet finite 
model. The absence of 1/s terms in Eqs. (45) and (46) is 
therefore a model artefact. 

We observe in addition that with attraction in the 
colour-octet channel; namely, for C > 0, (a^ — 1) is nec- 
essarily positive on the perturbative domain and a2 3 are 
negative. These results are also true in QCD (a^ > 1 
up to logarithmic corrections). It would be an excep- 
tional result if these statements were not also true on the 
nonperturbative domain. 

We now turn to the infrared domain and focus on s = 
but consider < m <C 1, in which case 



A{s = 0,to) 
B{s = 0,m) 



b^ + blm. 



(50) 
(51) 



Upon insertion of these expressions into the gap equation 
one obtains 



2 + 3C 
2 + C 



4 + 20C+25C^ 



6° = \/4 + 2C, bl 



(2+C)5/2 

1 

2 + C' 



(53) 



viz., results which show that in the neighbourhood of 

s = and with attraction in the colour-octet channel, 
A{s) decreases with increasing current-quark mass while 
B{s) increases. At this order the mass function is 



M(s = 0,to) 



S(0,to) 
A{0,m) 



with 



_b°o _ V2{2+Cr 



1 

,,1 - 1^ 



bl b°^a} 




1 

"0 



4-I-6C ' 
6 + 23C + 15C'^ 



(54) 

(55) 
(56) 



^0 W 2(2 + 3C)' 

For C > —(1/3) the mass function also increases with 
rising to. 

The infrared behaviour of the dressed-cjuark-gluon ver- 
tex follows immediately via Eqs. (36) - (38); using which 
one finds 



af(s = 0,TO) « afo + afi TO, z = 1,2,3, (57) 



with 



^1,0 



ai,i — 



"2' 
1 C 



^2,0 



^2,1 — 



=^3,0 



=^3,0 



V8 V2 + C' 

C(2-C)(2 + 3C) 
(2 + C)2 

2V2C (8 + 20C-f 2C2-9C3) 



(2 + C) 



7/2 



V8C 



V2 + C 
2C {5 + 6C) 
{2 + Cf 



(58) 
(59) 

(60) 

(61) 

(62) 

(63) 
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These algebraic formulae provide a clear indication of 
the effect on the dressed-quark-gluon vertex of attrac- 
tion in the projection of the quark-antiquark scattering 
kernel onto the colour-octet channel. Attraction causes 
ri{p,p) to be enhanced in the infrared cf. the bare ver- 
tex; and it drives i'2,3{p,p) to magnified, negative values. 
These results also signal that there are no surprises in the 
evolution into the infrared of the ultraviolet behaviour 
expressed in Eqs. (44) - (46): in each case attraction en- 
sures that a current-quark mass acts to reduce the vertex 
function in magnitude. 



B. Numerical results 

We have not hitherto specified a value for the param- 
eter C. Remember, a positive value ensures net attrac- 
tion in the colour-octet quark-antiquark scattering ker- 
nel and thereby models the dominance of the intrinsically 
non-Abelian £j contribution to the dressed-quark-gluon 
vertex over the Abelian-like term £^ . 

In choosing a value for C we elect to be guided by 
results from quenched lattice-QCD simulations of the 
dressed-quark propagator [28] and dressed-quark gluon 
vertex [45] . We focus on a current-quark mass common 
to both simulations; namely, 60 MeV, at which value the 
lattice dressed-quark propagator has [29]: .^qu(O) « 0.7, 
■Mqu(O) Ri 0.42. Then, so as to work with dimension- 



less quantities, we set me 



0.06/Mqu(0) and, using 



Eqs. (50) - (63), require a least-squares fit to^' 



A{meo) 


= 1.4, 


(64) 


ai(0,m6o) 


= 2.1, 


(65) 


-M{0,meof a^{0,m6o) 


= 7.1, 


(66) 


-M(0,m6o) 03(0, meo) 


= 1.0. 


(67) 


This procedure yields 






C = C = 0.51 




(68) 



with an average relative error f = 25% and standard 
deviation cr^ = 70%. It is plain that the fit is not perfect 
but, on the other hand, the model is simple. We note 
that for C = 0.6:, f = 21%, (t,. = 72%, while for C = 0.4: 
f = 31%, ar = 67%. If one omits Eq. (66) from the 
fitting requirements then C = 0.49 with f = 2.5% and 
ar = 63%. It is evident that competing requirements 
bound the amount of attraction allowed in the kernel. 
We can now illustrate the results for the dressed-quark 
propagator and dressed-quark-gluon vertex. 

In Figs. 1 3 we depict Aip"^), the scalar function 
which characterises the Dirac vector piece in the inverse 



[iv] There is a confusion of positive and negative signs in Ref. [45] 
concerning A2, A3, as defined therein. Our signs are correct. 
With the conventions expressed in Eq. (27): 4A2 = —0:2 and 
2 A3 = 03. 




FIG. 1: Current-quark- mass-dependence of the inverse of the 
drcsscd-quark wave function renormalisation. For all curves 
C = C = 0.51. Dash-dot line: m = 2; dotted line: m = meo; 
solid line: m = 0.015; dashed line: chiral limit, m = 0. All 
dimensioned quantities are measured in units of Q in Eq. (26). 
A fit to meson observables requires Q = 0.69 GeV and hence 
m = 0.015 corresponds to 10 MeV. 




FIG. 2: C-dependcncc of A{s). For all curves rn — 0.015. 
Solid line: C = C = 0.51; dash-dot-dot line: C = 1/4; dotted 
line: C — 0, which is equivalent to the rainbow truncation 
result first obtained in Ref. [33]; and dash-dot curve: C = 
— 1/8, for comparison with Ref. [32]. Dimensioned quantities 
are measured in units of G in Eq. (26). 



of the dressed-quark propagator. Figure 1 shows A{p'^) 
to evolve slowly with the current-quark mass when that 
mass is significantly smaller than the model's mass-scale. 
However, when the current-quark mass becomes com- 
mensurate with or exceeds that mass-scale, it acts to 
very effectively dampen this function's momentum de- 
pendence so that A{p'^) Ri 1. This is also true in QCD. 

Figure 2 exhibits the C-dependence of A{p'^). A com- 
parison of the C = —1/8, 0, 1/4 and 0.51 curves shows 



FIG. 3: Pointwise approach by solutions of the truncated gap 
equation to that obtained with the completely resummed ver- 
tex. Solid lino: complete solution; and in addition: dash- 
dasli-dot line - result for A{s) obtained with only the i — 0,1 
terms retained in Eq. (29) , which corresponds to the one- loop 
corrected vertex; short-dash line - two-loop-corrected vertex, 
i = 0,1, 2; long-dash line - three-loop-corrected vertex; and 
short-dash-dot line: four-loop corrected vertex. For all curves, 
C = C = 0.51 and m = 0.015. Dimensioned quantities are 
measured in units of Q in Eq. (26) . 



clearly that the presence of net attraction in the colour- 
octet quark-antiquark scattering kernel uniformly in- 
creases the magnitude of A{p'^) at all momenta. This 
effect is pronounced at infrared spacelikc momenta and 
particularly on the timelike domain, s < 0. In this fig- 
ure, as in those which follow, C = corresponds to the 
rainbow- ladder DSE tnmcation; i.e., the leading order 
term of the truncation scheme introduced in Ref. [31]. 

In Fig. 3 we illustrate the effect of adding corrections to 
the vertex, one at a time, in the presence of the amount 
of attraction suggested by lattice data. On the spacelike 
domain, .s > 0, the one-, two-, three- and four-loop cor- 
rected vertices yield a result for A{p^) that is little differ- 
ent from that produced by the completely resummed ver- 
tex (solid line). However, that is not true on the timelikc 
domain, whereupon confinement is expressed and hence 
it is unsurprising that nonperturbative effects become im- 
portant. In our model, as in the original version [33], 
confinement is realised via the absence of a particle-like 
singularity in the dressed-quark propagator [3] . The cusp 
displayed by A{p'^) in the timelike domain is one man- 
ifestation of this feature. While it moves toward s = 
with increasing C, it does not enter the spacelike domain. 
The figure makes clear that convergence to the solution 
obtained with the completely resummed vertex proceeds 
via two routes: one followed by solutions obtained with 
an odd number of loop corrections to the vertex; and 
another by those obtained using a vertex with an even 
number of loop corrections. This effect is not present 
with net repulsion in the colour-octet projection of the 
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FIG. 4: Upper panel - Current-quark- mass-dependence of the 
dressed-quark mass function. For all curves C — C — 0.51. 
Dotted line: m = meo; solid line: m = 0.015; dashed line: 
chiral limit, m = 0. Lower panel - C-dependence of this 
function. For all curves m = 0.015. Solid line: C = C = 0.51; 
dash-dot-dot line: C = 1/4; dotted line: C = 0, which is 
equivalent to the rainbow truncation result first obtained in 
Ref. [33]; and dash-dot curve: C — —1/8, for comparison with 
Ref. [32]. In addition, for C — 0.51: dash-dash-dot line - result 
for M(s) obtained with only the i = 0,1 terms retained in 
Eq. (29), which is the one- loop corrected vertex; and short- 
dash line - two-loop-corrected vertex. Dimensioned quantities 
are measured in units of Q in Eq. (26). 



quark-antiquark scattering kernel. 

We plot the dressed-quark mass function in Fig. 4. The 
upper panel illustrates the dependence of M(p^) on the 
current-quark-mass. The existence of a nontrivial solu- 
tion in the chiral limit is the realisation of DCSB, in our 
model and also in QCD. For current-quark masses less 
than the model's mass-scale, Q, the dynamically gener- 
ated mass determines the scale of observable quantities. 
However, for m}^ Q, this explicit chiral symmetry break- 
ing mass-scale overwhelms the scale generated dynami- 
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cally and enforces Mip^) w m. This is the behaviour of 
the 5-quark mass function in QCD [35]. 

In the lower panel of Fig. 4 we depict the C-dependence 
of M(p^). In discussing this figure we take the rainbow- 
ladder result, C = 0, illustrated with the dotted line, as 
our reference point. From a comparison of this curve with 
the C = —1/8 result (dash-dot hne) and the C = 0.51 
result (solid line) it is immediately apparent that ver- 
tex dressing driven by net attraction in the colour-octet 
quark-antiquark scattering kernel reduces the magnitude 
of the mass function at infrared momenta, a trend which 
is reversed for spacelike momenta s ^ Q}'"'^ This effect, 
which is correlated with the enhancement of A{p^), has 
an impact on the magnitude of the vacuum quark con- 
densate. 

For example, the mapping of Eq. (3) into our model is: 



87r2 



ds 8 



Z{s)M{s) 
s-HM(s)2 



(69) 



where sq is the spacelike point at which the model's mass 
function vanishes in the chiral limit, and we find 



-{m)c=c = (0.183 a)^ = (0.13 GeV)^ 



(70) 



with Q = 0.69 GeV. The rainbow-ladder result is 
-('?'?)c=o = e?V(207r2) = (0.172^)^ = (0.12 GeV)^ so 
that 



= 0.82. 



(71) 



This ratio drops to 0.50 when C = 1.0 is used to calculate 
the denominator. 

It is thus evident that in the presence of attraction 
in the colour-octet quark-antiquark scattering kernel and 
with a common mass-scale the condensate is significantly 
larger than that produced by a ladder vertex owing to an 
expansion of the domain upon which the dressed-quark 
mass function has nonzero support. 

It is natural to ask for the pattern of behaviour in the 
presence of repulsion. In this case Fig. 4 indicates that 
with C = — 1/8 the value of the mass function is enhanced 
at s = 0, a feature that is associated with the suppression 
of A{s) which is apparent Fig. 2. The magnitude of the 
mass function grows larger still with a further decrease in 
C and its domain of nonzero support expands. Therefore 
here, too, the condensate is larger than with the ladder 
vertex; e.g., 



0.92, 



(72) 



[v] This pattern of behaviour is familiar from explorations [20, 22] of 
the effect in the gap equation of vertex Ansdtze [13, 46]; i.e., mod- 
els for the vertex whose diagrammatic content is unknown but 
which exhibit properties in common with our calculated C > 
result. 



and the ratio drops to 0.49 when C = —3/8 is used to 
evaluate the denominator. 

The implication of these results is that in general, with 
a given mass-scale and a common model dressed-gluon 
interaction, studies employing the rainbow-ladder trun- 
cation will materially underestimate the magnitude of 
DCSB relative to those that employ a well-constrained 
dressed-quark-gluon vertex. Naturally, in practical phe- 
nomenology, alterations of the mass-scale can compen- 
sate for this [20]. 

Figure 4 also shows that the pointwise convergence of 
M(s) to the solution obtained with the fully resummed 
vertex, by solutions obtained with truncated vertices, fol- 
lows a pattern similar to that exhibited by A{s). 

In the upper panel of Fig. 5 we illustrate the current- 
quark-mass-dependence of the scalar function associated 
with 7^ in the dressed-quark-gluon vertex. In QCD this 
is the only function that displays an ultraviolet diver- 
gence at one- loop in perturbation theory. As with A{p^), 
c^iip'^) evolves slowly with the current-quark mass but 
again, when the current-quark mass becomes commen- 
surate with or exceeds the theory's mass-scale, it acts to 
very effectively dampen the momentum dependence of 
this function so that a;i(p^) « 1. This effect is appar- 
ent in the rainbow vertex model employed in Ref. [29] to 
describe quenched-QCD lattice data. 

The lower panel of Fig. 5 portrays the C-dependence 
of ai{p'^). It is here particularly useful to employ the 
rainbow- ladder result, C = 0, illustrated with the dot- 
ted line, as our reference point, because this makes the 
contrast between the effect of attraction and repulsion in 
the colour-octet quark-antiquark scattering kernel abun- 
dantly clear. Attraction uniformly increases the magni- 
tude of ai(p2), while the opposite outcome is produced 
by omitting the effect of the three-gluon-vertex in the 
DSE for the dressed-quark-gluon vertex. 

The upper panel of Fig. 6 illustrates the current-quark- 
mass-dependence of the scalar function modulating the 
subleading Dirac vector component of the dressed-quark- 
gluon vertex. For light-quarks this function is material in 
size and negative in magnitude. However, its importance, 
too, diminishes with increasing current-quark mass. In 
the lower panel we show the C-dependence of a'^ip'^). 
The qualitative features of the completely resummed re- 
sult for a'i{p'^) are also manifest here. However, for this 
component of the vertex, which is purely dynamical in 
origin, there is a marked difference at timelike momenta 
between the result obtained with an odd number of loop 
corrections in the vertex and that obtained with an even 
number. Here the dual pathway pointwise convergence 
to the fully resummed result is conspicuous. 

It is notable that the size of our complete result for 
a'^ip'^) is an order of magnitude smaller than that re- 
ported in Ref. [45]. This is an isolated case, however. 
The calculated magnitudes of the other functions in the 
dressed-quark propagator and dressed-quark-gluon ver- 
tex are commensurate with those obtained in quenched- 
QCD lattice simulations. We remark in addition that 



10 




I Q5 I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

-4 -2 2 4 



1.2 - 



1 - 




I 1 I 1 I 1 I 1 I 1 I 1 I 

-4 -2 2 4 6 8 

s 



FIG. 5: Upper panel - Current-quark-mass-dependence of 
Qi (s); viz., the leading component in the dressed-quark-gluon 
vertex. For all curves C = C = 0.51. Dash-dot line: m = 2; 
dotted lino: m = mfio; solid line: m = 0.015; dashed line: 
chiral limit, rn = 0. Lower panel C-dcpendencc of this func- 
tion. For all curves rn = 0.015. Solid line: C — C = 0.51; 
dash-dot-dot line: C = 1/4; dotted line: C = 0, namely, the 
bare vertex, which defines the rainbow truncation; and dash- 
dot curve: C = —1/8, for comparison with Ref. [32]. In addi- 
tion, for C = 0.51: dash-dash-dot line - qi(s) obtained with 
only the t = 0, 1 terms retained in Eq. (29), which is the one- 
loop corrected vertex; and short-dash line - two-loop-corrected 
vertex. Dimensioned quantities are measured in units of Q in 
Eq. (26). 
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FIG. 6: Upper panel - Current-quark-mass-dependence of 
a^{s) in Eq. (27). For all curves C = C = 0.51. Dash-dot line: 
m = 2; dotted line: m = meo; solid line: m = 0.015; dashed 
line: chiral limit, m = 0. Lower panel - C-dcpendcncc of this 
function. For all curves rn — 0.015. Solid line: C — C = 0.51; 
dash-dot-dot line: C = 1/4; and dash-dot curve: C = —1/8, 
for comparison with Ref. [32]. Moreover, for C = 0.51: dash- 
dash-dot line - one- loop result for a^is), namely, the result 
obtained with only the i = 0,1 terms retained in Eq. (29); 
short-dash line - two- loop-corrected result; long-dash line - 
three-loop-corrected; and short-dash-dot line: four-loop cor- 
rected. NB. There is no C = curve because a2 («) = in the 
bare vertex. Dimensioned quantities are measured in units of 
g in Eq. (26). 



the lattice result is an order of magnitude larger than 
that obtained with a commonly used vertex Ansatz [46] . 
This curious discrepancy deserves study in more realistic 
models. 

In Fig. 7 we display the result for what might be called 

the scalar part of the drcsscd-quark-gluon vertex; viz., 
a'^{p^). This appellation comes from the feature, ap- 
parent in Eq. (46), that a'^ip^) is the piece of the ver- 
tex whose ultraviolet behaviour is most sensitive to the 
current-quark mass. The figure demonstrates, too, that 



at infrared momenta a^{p^) is materially affected by the 
scale of DCSB, Eq. (70): at s = the deviation from 
its value in the rainbow truncation is approximately four 
times that exhibited by aj(p^). Hence, this term can be 
important at infrared and intermediate momenta. How- 
ever, for light quarks, owing to its sensitivity to the 
current-quark-mass, it rapidly becomes unimportant at 
ultraviolet momenta. The dual pathway pointwise con- 
vergence to the fully rcsummed result in the timelike re- 
gion is also prominent here. 
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FIG. 7: Upper panel - Currcnt-quark-mass-dependence of 
af (s) in Eq. (27). For all curves C = C = 0.51. Dash-dot line: 
m = 2; dotted line: m = meo; solid line: m = 0.015; dashed 
line: chiral limit, m = 0. Lower panel - C-dcpcndcncc of this 
function. For all curves rn = 0.015. Solid line: C — C = 0.51; 
dash-dot-dot line: C — 1/4; and dash-dot curve: C = —1/8, 
for comparison with Ref. [32]. Moreover, for C = 0.51: dash- 
dash-dot line - one- loop result for 0.3(3), namely, the result 
obtained with only the i = 0,1 terms retained in Eq. (29); 
short-dash line - two- loop-corrected result; long-dash line - 
thrcc-loop-corrcctcd; and short-dash-dot line: four-loop cor- 
rected. NB. There is no C = curve because (s) = in the 
bare vertex. Dimensioned quantities are measured in units of 
g in Eq. (26). 



Finally, since they are absent in rainbow truncation, 
it is illuminating to unfold the different roles played by 
Q!2 (s) and ag (s) in determining the behaviour of the self- 
consistent solution of the gap equation. Some of these ef- 
fects are made evident in Fig. 8 and all can readily be un- 
derstood because Eqs. (40) and (41) signal much of what 
to expect. The key observation is that a'^{s) alone is 
the source of all coupling between Eqs. (40) and Eq. (41) 
that is not already present in rainbow-ladder truncation: 
a'^ B appears in the equation for A[s) and a'^ s A appears 
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FIG. 8: Upper panel - Impact of ol^^s) and ok,(s) on A{s). For 
C = C, Eq. (68), solid line: result obtained with both terms 
present in the vertex; dashed-line: 02(5) omitted; dash-dot 
line: 03(3) omitted. The dotted line is the result obtained 
with both terms present in the vertex but C = —1/8. Lower 
panel - Impact of 0:2(5) and 0:3(5) on M(s). The legend is 
identical. In all cases m = 0.015 and dimensioned quantities 
are measured in units of Q in Eq. (26) . 



in the equation for B{s). The action of is merely to 
modify the rainbow-ladder coupling strengths. 

Turning to detail, each equation involves the combi- 
nation 50:2(3), which indicates that the impact of 0:2(5) 
on both A(s) and B(s) is negligible at infrared momenta 
but possibly material on the intermediate spacelikc do- 
main. This suppression is also felt by 03(5) in Eq. (41) 
but not in Eq. (40). Hence, omitting 03(5) will affect 
A{s) at infrared momenta but not B{s). These expecta- 
tions are easily substantiated by repeating the analysis 
that gave Eqs. (52), (53). When of = is imposed in 
the gap equation one finds Oq = 2; namely, the same re- 
sult that is obtained in the rainbow truncation, but no 
change in 6q. For C > this naturally leads to an increase 
in M(0) = b^/a^. On the other hand, fixing a^{s) = 
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produces no change in aj] or b^. These resuhs are appar- 
ent in the figure. (NB. It is a numerical accident that 
the dash-dot and dotted curves in the lower panel almost 
coincide at s = 0.) 

It will readily be appreciated that neither a2{s) nor 
03(5) can affect the deep ultraviolet behaviour of the gap 
equation's solution, Eqs. (42), (43), (47), because they 
vanish too rapidly as 1/s ^ 0. This is plain in Fig. 8. 

The effects at intermediate spacelike momenta can also 
be inferred by inspection of Eqs. (40), (41). Both a2{s) 
and a^{s) are negative and hence serve to magnify A{s) 
with respect to the value obtained using a bare vertex. 
However, they compete in Eq. (41), with a2{s) work- 
ing to diminish -B(.s) and 0:3 (.s) acting to amplify it. 
In the absence of a'^{s), therefore, one should expect 
M(s) = B{s)/A{s) to be suppressed at intermediate mo- 
menta because the numerator is reduced and the denom- 
inator increased, and consequently a condensate much 
reduced in magnitude. Omitting 02(5) should yield the 
opposite effect. This is precisely the outcome of our nu- 
merical studies: 



c=o 



c=c 



= 1.97; and 



c=o 



c=c 



0.40 . 



(73) 

At timelike momenta the omission of 0:2(5) has little 
qualitative impact on A{s) and M(s) because, as we have 
seen, 0:3(5) is always the dominant term at infrared mo- 
menta and hence provides the primary seed for evolution 
to s < 0. In the absence of 03 (s), on the other hand, the 
momentum dependence of A{s) and M{s) must resem- 
ble that obtained in rainbow truncation because 502(5) 
has only a small domain of material support. Figure 8 
makes this clear. (NB. The near pointwise agreement of 
the dash-dot and dotted curves in the lower panel is a 
numerical accident.) Importantly, the confining charac- 
ter of the dressed-quark propagator is preserved in both 
cases. Indeed, omitting both 03(5) and 03(5) simultane- 
ously does not qualitatively alter this feature of the model 
prescribed by Eq. (26), although the omission does affect 
aspects of its expression. That, however, is unsurprising. 
It is evident, for example, in a comparison of Refs. [21, 22] 
with Ref. [33]. 



IV. BETHE-SALPETER EQUATION 

The renormalised homogeneous Bethe-Salpeter equa- 
tion (BSE) for the quark-antiquark channel denoted by 
M can be compactly expressed as 

.A 

[TM{k;P)]EF= [K{k,q;P)]'i^[xM{q;P)]GH (74) 

Jq 

where: k is the relative momentum of the quark- 
antiquark pair and P is their total momentum; E, . . . ,H 
represent colour, fiavour and spinor indices; and 



with rM{q',P) the meson's Bethc-Salpcter amplitude. 
In Eq. (74), K is the fully-amputated dressed-quark- 
antiquark scattering kernel. 



A. Vertex consistent kernel 

The preservation of Ward-Takahashi identities in those 
channels related to strong interaction observables re- 
quires a conspiracy between the dressed-quark-gluon 
vertex and the Bethe-Salpeter kernel [31, 47]. The 
manner in which these constraints are realised for ver- 
tices of the class considered herein was made explicit 
in Ref. [32]. In that systematic and nonperturbative 
truncation scheme the rainbow gap equation and lad- 
der Bethe-Salpeter equation represent the lowest-order 
Ward-Takahashi identity preserving pair. Beyond this, 
each additional term in the vertex generates a unique 
collection of terms in the Bethe-Salpeter kernel, a subset 
of which are always nonplanar. 

For any dressed-quark-gluon vertex appearing in the 
gap equation, which can be represented expressly by an 
enumerable series of contributions, the Bethe-Salpeter 
kernel that guarantees the validity of all Ward-Takahashi 
identities is realised in 



TM{k;P)= [ V^^{k-q)r^^ 

X \xM{q; P) r,(g_, fc_) + S{q+) Alj,{q, k; P) 



where 



(76) 



A^,(g,fc;P) = ^AX(9,fc;P), (77) 

n=0 



with 



XM{k; P) = S{k+) TM{k; P) S{k-) , 



(75) 



AS(^, fc; P) = fvp,{l - q) l\ XMiq; P) 

X (g_, g_ + k-i) S{q- + k - £) / V 

.A 

+ / Vp4k-q)l'''ypS{q++i-k) 

Jq 

X '"r^,n-l (9+ +^-k,q+) XMiq; P) I'la 

+ I Vp,{£-q')HS{q'+) 

Jq' 

X A»'"-i(g', q' + k-£;P) 5(g'_ + k - £) Z^. 

(78) 

This last equation is a recursion relation, which is to be 
solved subject to the initial condition A^°^ = 0. 

It is apparent that the kernel of Eq. (76) is completely 
known once the vertex to be used in the gap equation is 
fixed and that equation's solution is obtained. 



13 



The Bcthc-Salpeter amplitude for any meson can be 
written in the form 



Nm 



TM{k; P)=IcYl ^'(^; ^) fMik\k-P; P^) =: [G] , 

(79) 

whore Q^{k;P) are those independent Dirac matrices re- 
quired to span the space containing the meson under con- 
sideration. It then follows upon substitution of this for- 
mula that Eq. (78) can be written compactly as 



.a;« 



4;n-lJ 



I f M ■ 



(80) 



This expression states that A^" can primarily be con- 
sidered as a matrix operating in the space spanned by 
the independent components of the Bethe-Salpeter am- 
plitude, with its Dirac and Lorentz structure projected 
via the contractions in the BSE. The first term (/Cm) in 
Eq. (80) represents the contribution from the first two in- 
tegrals on the right-hand-side (r.h.s.) of Eq. (78). This 
is the driving term in the recursion relation. The second 
term {£) represents the last integral, which enacts the 
recursion. 



B. Solutions of the vertex-consistent meson 
Bethe-Salpeter equation 

A discussion of the general features of the vertex- 
consistent BSE and its sohition is given in Ref. [32]. 
Herein we therefore focus immediately on the model un- 
der consideration. 



1. n -meson 

With the model of the dressed-gluon interaction in 
Eq. (26) the relative momentum between a meson's con- 
stituents must vanish. It follows that the general form of 
the Bethe-Salpeter amplitude for a pseudoscalar meson 
of equal- mass constituents is 



r.(P) = 75 ifi{P^) + i-Pf2{P^) 



(81) 



where P, = 1, is the direction- vector associated with 
P. 

To obtain the vertex-consistent BSE one must first de- 
termine A°',}. That is obtained by substituting Eq. (81) 
into the r.h.s. of Eq. (78). Only the first two integrals 
contribute because of the initial condition and they are 
actually algebraic expressions when Eq. (26) is used. This 
gives /CJry Eq. (80). Explicit calculation shows this to 
be identically zero, and hence A"'^ = 0. Furthermore 
because this is the driving term in the recursion relation 
then 



The result is not accidental. It is fundamentally a con- 
sequence of the fact that A^^(£, A:;P) = /" A^^(f , /e; P) 
is an axial- vector vertex which transforms under charge 
conjugation according to: 



A,^(-fc,-^;P)' = -A,^(^,fc;P), 



(83) 



where (•)* denotes matrix transpose, and Eq. (26) en- 
forces k = Q = I m. the BSE. It is not a general feature 
of the vertex-consistent Bethe-Salpeter kernel. 

One thus arrives at a particularly simple vertex- 
consistent BSE for the pion (Q=P/2): 

r.(P) = - 7^ S{Q) r,(P) S{-Q) r^(-Q) . (84) 



Consider the matrices 



'Pi = -t75, P2 = -,1-Plh- 



which satisfy 



f:{P)=tr^ViT^{P). 



(85) 



(86) 



They may be used to rewrite Eq. (84) in the form 

f^{p) = nAP')f^{p), (87) 

wherein H„{P'^) is a 2 x 2 matrix 

'^AP^)i3 = 

-^troP, 7m s{Q) r.(p) si-Q) r^(-g) . 

(88) 

Equation (87) is a matrix eigenvalue problem in which 
the kernel W is a function of P^. This equation has a 
nontrivial solution if, and only if, at some 

det[n^{P')-l]\p,^^,^^ = 0. (89) 

The value of M for which this characteristic equation is 
satisfied is the bound state's mass. NB. In the absence of 
a solution there is no bound state in the channel under 
consideration. 

We have solved Eq. (89) for the pion and the results 
are presented in Table I. The fact that the vertex- 
consistent Bethe-Salpeter kernel ensures the preservation 
of the axial-vector Ward-Takahashi identity, and hence 
guarantees the pion is a Goldstone boson in the chiral 
limit, is abundantly clear: irrespective of the value of C 
and the order of the truncation, = for m = 0. Away 
from this symmetry-constrained point the results indi- 
cate that, with net attraction in the colour-octet quark- 
antiquark scattering kernel, the rainbow-ladder (lowest 
order) truncation overestimates the mass; i.e., it yields a 
value greater than that obtained with the fully resummed 
vertex (n = oo). Moreover, the approach to the exact re- 
sult for the mass is not monotonic. On the other hand, 
given two truncations for which solutions exist, charac- 
terised by ni- and n2-loop insertions, respectively, then 



n2 > ni ; (90) 



Ae,(A;,fc;P)^0. 



(82) 



viz., correcting the vertex improves the accuracy of the 
mass estimate. 
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TABLE I: Calculated tt and p meson masses, in GeV, quoted 
with g = 0.69 GeV, in which case m = 0.015 5 = lOMeV. 
(In the notation of Ref. [31], this value of G corresponds to 
rj = 1.39 GoV.) n is the number of loops retained in dress- 
ing the quark-gluon vertex, see Eq. (29), and hence the order 
of the vertex-consistent Bethe-Salpeter kernel. NB. n = 
corresponds to the rainbow-ladder truncation, in which case 
rup = V^g, and that is why this column's results are inde- 
pendent of C. 
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TT, 
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= 0.01 


0.149 
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0.140 


0.138 






P, 
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= 


0.982 




0.828 


0.754 






P, 


m 


= 0.01 


0.997 




0.844 


0.770 



where the projection matrices are 

^1^ = ^ 7 • e\P) ,r^ = ^ <T^. e^(P) Pu . (96) 



The remaining term in Eq. (94) is generated by A^^, 
which is nonzero in this case, and produces planar and 
nonplanar contributions. To be explicit 

A%ie\P) = l-Ap,Ae\P), (97) 

where the Dirac structure is completely expressed 
through 

Ap4e\ P) = /3i(P)e^ + (32{P) eth- P 

+ l3s{P) ij ■e^P^+ l3i{P) <7„^ e^(P) 7, 

+ /35 (P) i<Joc0 el (P) P/3 P. + /36 (P) 7 • eS. • (98) 

The manner by which the vector /3(P) = 
column[/3i,/32,/33,/34,/35,/36] is determined was de- 
scribed in Ref. [32] and in our case that procedure 
yields 



1 



1 - L,{P) 



Gp{P) 



(99) 



with 



2. p-meson 

In our algebraic model the complete form of the Bethe- 
Salpeter amplitude for a vector meson is 

r^(P) = 7 • e\P) /f (P2) + a,. e^(P) P. f^{P') . (91) 

This expression, which has only two independent func- 
tions, is miich simpler than that allowed by a more realis- 
tic interaction, wherein there are eight terms. Neverthe- 
less, Eq. (91) retains the amplitudes that are found to be 

dominant in more sophisticated studies [48]. In Eq. (91), 
{e^(P); A — — 1,0,+1} is the polarisation four- vector: 



P-e^{P) = 0,yX; e^{P)-e^'{P) = 6 



rAA' 



(92) 



The construction of the vertex-consistent BSE for the 
class of vertices under consideration herein is fully de- 
scribed in Ref. [32] . With the model for the dressed-gluon 
interaction in Eq. (26) , that procedure yields an equation 
that can be written 



/,(P)=H,(P^)/,(P), 
wherein Wp(P^) is a 2 x 2 matrix: 

Hp{P^) =HpT{P^) +HpA{P^ 



(93) 



(94) 



In Eq. (94), HpriP^) is analogous to the sole term that 
remains for the pion; viz., 

Wpr(P^)ij = 

- -^truV,^ 7^ S{Q) r^(P) S{-Q) r^(-g) , (95) 



-^^'LpiP) 



2A 
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-2A 
2P2 

A 








-2QAB 






2A 


-2QAB 






(100) 



in which the argument of each function is Q^, and 
Gp{P) = 

2[fP{P)B{Q^) + fpQA{Q^^^ 



-4Ca?(Q) 
A(Q) 



/f(P)Q^(Q2)^/P(P)ij(g2) 

-f,{P)QA{Q^) + fmB{Q^) 
-f?iP)QA{Q^)+f^{P)B{Q^) 







(101) 



It is correct that in this study G p{P) depends only on 
and, moreover, the zeros in its last two rows entail 

/35(P)-0 = /36(P). 

The last term in Eq. (94) is thus 



W,a(P') = 



_5_ 



ptvj,Vt^^S{Q)Kp^{P) (102) 



and hence the p-meson's mass, Mp, is obtained as the 
solution of 



= det [W^(P2) - /] 



p2 + M2=0 



(103) 



det[W^r(P')+HpA(P')-/] 
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We have solved Eq. (103) and present the results in 
Table I. One sees that with increasing net attraction in 
the colour-octet projection of the quark- antiquark scat- 
tering kernel the amount by which the rainbow-ladder 
truncation overestimates the exact mass also increases: 
with the amount of attraction suggested by lattice data 
the n = mass is 27% too large. A related observation 
is that the bound state's mass decreases as the amount 
of attraction between its constituents increases. This is 
unsurprising. Furthermore, with increasing attraction, 
even though the fully resummed vertex and consistent 
kernel always yield a solution, there is no guarantee that 
a given truncated system will support a bound state: the 
one- loop corrected vertex and consistent kernel (n = 1) 
do not have sufficient binding to support a /O-meson. This 
is overcome at the next order of truncation, which yields 
a mass 9.7% too large. Given the results for the dressed- 
quark propagator and dressed-quark-gluon vertex illus- 
trated in Figs. 1-7; namely, the dual pathways of point- 
wise convergence to the true solutions, the difference be- 
tween odd- and even-n truncations in our model is under- 
standable. Nevertheless, the primary observation; viz., 
a given beyond-rainbow-ladder truncation may not sup- 
port a boimd state even though one is present in the so- 
lution of the complete and consistent system, provides a 
valuable and salutary tip for model building and hadron 
phenomenology. 

Finally, as has often been observed, and independent 
of the truncation, bound state solutions of gap-equation- 
consistent BSEs always yield the full amount of n-p mass 
splitting, even in the chiral limit. This splitting is plainly 
driven by the DCSB mechanism and therefore its true un- 
derstanding requires a veracious realisation of that phe- 
nomenon. 



3. Dependence on the current-quark mass 

In connection with this last observation it is relevant 
to explore the evolution with current-quark mass of the 
pseudoscalar and vector meson masses and, consequently, 
of the difference between them. The results for pseu- 
doscalar mesons should be interpreted with the following 
caveat in mind. In constructing the vertex and kernel 
we omitted contributions from gluon vacuum polarisation 
diagrams. These contribute only to flavour singlet me- 
son channels. Hence, for light-quarks in the pseudoscalar 
channel, wherewith such effects may be important [49— 
51], our results should be understood to apply only to 
flavour nonsinglets. In principle, the same is true for 
light vector mesons. However, experimentally, the lo and 
(j) mesons are almost ideally mixed; i.e., the a; exhibits no 
ss content whereas the cj) is composed almost entirely of 
this combination. We therefore assume that the vacuum 
polarisation diagrams we have omitted are immaterial 
in the study of vector mesons. (NB. It is an artefact 
of Eq. (26) that this model supports neither scalar nor 
axial- vector meson bound states [32, 33].) 



TABLE II: Current-quark masses required to reproduce the 

experimental masses of the vector mesons. The values of m,,^,, 
m,,j_ arc predictions. Experimentally [52], m,,^ = 2.9797 ± 
0.00015 and m^, = 9.30 ± 0.03. NB. O^j is a fictitious pseu- 
doscalar meson composed of unlikc-flavour quarks with mass 
nis, which is included for comparison with other nonpertur- 
bative studies. All masses are listed in GeV. 





= 0.01 


rus = 


= 0.166 


rric 


= 1.33 




= 4.62 


nip 


= 0.77 




= 1.02 




= 3.10 


mr(is) 


= 9.46 




= 0.14 




= 0.63 




= 2.97 




= 9.42 



We fix the model's current-quark masses via a fit to the 
masses of the vector mesons, and the results of this exer- 
cise are presented in Table II. The model we're employing 
is ultraviolet finite and hence our current-quark masses 
cannot be directly compared with any current-quark 
mass-scale in QCD. Nevertheless, the values are quantita- 
tively consistent with the pattern of flavour-dependence 
in the explicit chiral symmetry breaking mass-scales of 
QCD. 

Our calculated results for the current-quark mass- 
dependence of pseudoscalar and vector meson masses are 
presented in Fig. 9. In the neighbourhood of the chiral 
limit the vector meson mass is approximately indepen- 
dent of the current-quark mass whereas the pseudoscalar 
meson mass increases rapidly, according to (in GeV) 

ml- w 1.33 m m<^g, (104) 

thereby reproducing the pattern predicted by QCD [34] . 

With the model's value of the vacuum quark conden- 
sate, Eq. (70), this result allows one to infer the chiral- 
hmit value of /° = 0.056 GeV via the Gell-Mann-Oakes- 
Renner relation. It is an artefact of the model that the 
relative-momentum-dependence of Bethe-Salpeter ampli- 
tudes is described by (5^(p) and so a direct calculation of 
this quantity is not realistic. The value is low, as that 
of the condensate is low, because the model is ultraviolet 
finite. In QCD the condensate and leptonic decay con- 
stant are much influenced by the high-momentum tails 
of the dressed-quark propagator and Bethe-Salpeter am- 
plitudes. 

The curvature in the pseudoscalar trajectory persists 
over a signiflcant domain of current-quark mass. For ex- 
ample, consider two pseudoscalar mesons, one composed 
of unlike-flavour quarks each with mass 2ms and another 
composed of such quarks with mass mg. In this case 

ml_ 

^ = 2.4, (105) 

a result which indicates that the nonlinear evolution ex- 
hibited in Eq. (104) is still very much in evidence for 
current-quark masses as large as twice that of the s- 
quark. With this result our simple model reproduces a 
feature of more sophisticated DSE studies [53-55] and a 
numerical simulation of quenched lattice-QCD [57]. 
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FIG. 9: Evolution of pseudoscalar and vector qq meson masses 
with the current-quark mass. Solid lino: pseudoscalar me- 
son trajectory obtained with C = C — 0.51, Eq. (68), using 
the completely resummed dressed-quark-gluon vertex in the 
gap equation and the vertex-consistent Bethe-Salpeter kernel; 
short-dash line: this trajectory calculated in rainbow-ladder 
truncation. Long-dash line: vector meson trajectory obtained 
with C using the completely resummed vertex and the con- 
sistent Bethe-Salpotor kernel; dash-dot line: rainbow-ladder 
truncation result for this trajectory. The dotted vertical hues 
mark the current-quark masses in Table 11. 

The mode of behaviour just described is overwhelmed 
when the current-quark mass becomes large: Q. In 

this limit the vector and pseudoscalar mesons become de- 
generate, with the mass of the ground state pseudoscalar 
meson rising monotonically to meet that of the vector 
meson. In our model 

^ = 1.04, (106) 

with a splitting of 130 MeV, and this splitting drops to 

just 40 MeV at mj; viz., only 5% of its value in the chiral 
limit. In addition to the calculated value, the general pat- 
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FIG. 10: Evolution with current-quark mass of the differ- 
ence between the squarcd-niasses of vector and pseudoscalar 
mesons. These results were obtained with C = 0.51, Eq. (68), 
using the completely resummed dressed-quark-gluon vertex 
in the gap equation and the vertex-consistent Bethe-Salpeter 
kernel. The dotted vertical lines mark the current-quark 
masses in Table 11. 



tern of our results argues for the mass of the pseudoscalar 
partner of the T{1S) to lie above 9.4 GeV. Indeed, we ex- 
pect the mass splitting to be much less than m j/^ — m^^ , 
not more. (See also; e.g., Ref. [56].) 

In Fig. 10 we depict the evolution with current-quark 
mass of the difference between the squared-masses of vec- 
tor and pseudoscalar mesons. On a material domain 
of current-quark masses this difference is approximately 
constant: to^_ — m^- « 0.56 GeV^, an outcome consis- 
tent with experiment that is not reproduced in numerical 
simulations of quenched lattice-QCD [57]. The difference 
is maximal in the vicinity of rUc, a result which reempha- 
sises that heavy-quark effective theory is not an appro- 
priate tool for the study of c-quarks [35] . 

Figure 11 is instructive. It shows that with growing 
current-quark mass the rainbow-ladder truncation pro- 
vides an increasingly accurate estimate of the ground 
state vector meson mass. At the s-quark mass the rel- 
ative error is 20% but that has fallen to < 4% at the 
c-quark mass. Similar statements are true in the valid 
pseudoscalar channels. In fact, in this case the agree- 
ment between the truncated and exact results is always 
better; e.g., the absolute difference reaches its peak of 
« 60 MeV at TO ~ 4 rris whereat the relative error is only 
3%. This behaviour in the pseudoscalar channel is funda- 
mentally because of Goldstone's theorem, which requires 
that all legitimate truncations preserve the axial-vector 
Ward-Takahashi identity and hence give a massless pseu- 
doscalar meson in the chiral limit. It is practically use- 
ful, too, because it indicates that the parameters of a 
model meant to be employed in a rainbow-ladder trun- 
cation study of hadron observables may reliably be fixed 
by fitting to the values of quantities calculated in the 
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current-quark-mass (GeV) 

FIG. 11: Evolution with current-quark mass of the differ- 
ence between the meson mass calculated in the rainbow-ladder 
truncation and the exact value; namely, that obtained using 
the completely resummed dressed-quark-gluon vertex in the 
gap equation and the vertex-consistent Bethe-Salpeter ker- 
nel. The upper panel depicts the absolute error and the lower 
panel, the relative error. Solid lines: vector meson trajecto- 
ries; and dashed-lines; pseudoscalar meson trajectories. The 
results were obtained with C = 0.51, Eq. (68) and the dotted 
vertical lines mark the current-quark masses in Table II. 



C. Bethe-Salpeter equation for diquarks 

Colour antitriplet quark-quark correlations (diquarks) 
have long been a focus of attempts to understand baryon 
structure [58, 59]. It was quickly realised that both 
Lorentz scalar and axial- vector diquarks, at least, are im- 
portant for baryon spectroscopy [60] and, from a consid- 
eration of baryon magnetic moments, that diquark cor- 
relations are not pointlike [61]. An appreciation of their 
importance has grown and a modern picture of diquark 
correlations in baryons is realised through their role in a 
Poincare covariant Faddeev equation [4, 62-64] . 

In rainbow-ladder truncation colour-antitriplct di- 
quarks are true bound states [65, 66]. (NB. The other ad- 
missible channel; viz., colour sextet, is never bound since 
even single gluon exchange is repulsive therein.) In spite 
of this, the addition of C^ik^p), Eq. (17), to the quark- 
gluon vertex, along with the three terms it alone gener- 
ates in the colour-antitriplet quark-quark scattering ker- 
nel, overwhelms the attraction produced by single gluon 
exchange and eliminates diquarks from the spectrum 

[31, 67]. The repulsive effect owing to £2 (^)P) is consum- 
mated when the series it generates is fully resummed [32]. 
With that vertex, and the consistent colour-antitriplet 
quark-quark scattering kernel, the characteristic polyno- 
mial obtained from the Bethe-Salpeter equation exhibits 
a pole, which is the antithesis of the zero associated with 
a bound state. 

Nevertheless, herein we have introduced a new element 
into the consideration of diquark correlations; namely, 
our model for the colour-octet quark-antiquark scatter- 
ing kernel exhibits attraction, which is a property not 
possessed by £2 {k,p)- The dressed-quarks that appear in 
the diquark Bethe-Salpeter equation are described by the 
gap equation we have already elucidated. However, the 
manner in which diagrams are combined and resummed 
in the vertex-consistent colour-antitriplet diquark Bethe- 
Salpeter kernel is different from that which maintains 
for colour singlet mesons. Fortunately, the modifications 
necessary in the class of Ansdtze containing our model 
were elucidated in Ref . [32] and we need only adapt them 
to our particular case. 

In consequence we arrive at the following Bethe- 
Salpeter equation: 



neighbourhood of the chiral limit. 

The general observation suggested by this figure is 
that with increasing current-quark mass the contribu- 
tions from nonplanar diagrams and vertex corrections 
are suppressed in both the gap and Bethe-Salpeter equa- 
tions. Naturally, they must still be included in precision 
spectroscopic calculations. It will be interesting to re- 
analyse this evolution in a generalisation of our study to 
mesons composed of constituents with different current- 
quark masses, and thereby extend and complement the 
limited such trajectories in Refs. [53, 54]. 



= -^7,S{Q)T^^{P)S{-Q)r''^{-Q) 
+^7m5(Q)T^^(P) (107) 

wherein: T^q{P) = Tqg{P)C'^, C = 7274 is the charge 
conjugation matrix, with Tqq{P) the Bethe-Salpeter am- 
plitude describing diquark correlations; and S and are 
precisely the elements used in the tt- and p-meson equa- 
tions. T^{P) is new, however. It is a variant for diquarks 
of in Eq. (77) , as we now show. 

The Bethe-Salpeter amplitudes for scalar and axial 
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vector diquarks can be written: 



that 



C1+ 



75 [A" {P')+1-Pf^ {P')\ , (108) 

(P) = 7-e^P)/r(P')+a^.e^(P)P./r(P'). 

(109) 



These ampHtudes are identical in structure to those of the 

TT- and /9-mcsons, Eqs. (81) and (91), and that is why we 
work with instead of r^^. Consequently, the Bethe- 
Salpeter equations assume the form in Eq. (93) 



wherein Ti,qq{P^) is a 2 x 2 matrix: 



(110) 



Wgg(P') = W,,r(P') + WggA(P') • (HI) 

The ?iqqr(P^) term is obtained in a straightforward man- 
ner from the first line of Eq. (107) by following the pat- 
tern of Eq. (95). 

The second term employs the mapping of into the 
present case, which is 



(112) 



n=0 



with 



—A"'" (P) 
6C '^"•'^ ' 



= i'ip><'kx^,ip)^in-ihQ) rr si-Q) {i'Tj, 

- 1\ S{Q) rr^„_,(g) x'^^XqqiP) {I'fia 

- l\ S{Q) A^-"-i(P) S{-Q) {I^Yj, , (113) 
where the superscript "t" denotes matrix transpose, 

xZ{P)=S{Q)r^qiP)S{-Q) (114) 

and {A^ = f, = l^, A^ = l"^} are three antisymmetric 
matrices that represent the colour structure of an an- 
titriplet diquark. For clarity we have already made use 
of Eq. (26). 

In general 

A«^^(P) = AUP) + A^,(P) A^ , (115) 



1 

8C 



7'L;n \k 



:= -^rA 
8C 



Du 



^[2AL-A^JA^ 

K 1 

(-11,11-1/ p\ r<'^ 



HP) 



1 <j 



5 

12 



(116) 



-L nrR-n \ k 



1 A a;n /;a\t 



1 

12C 
1 

12" ^99!- 

5 



A?, + 2A^JA^ 
5 

12 



^Gl'2-\P) + Y^Gll-HP) 

^I'ql ^ (P) ~ Y2 "^99" ^ ' 



wherein 



(P) 



qgu 

G%l{P) 



(117) 

lpXqg{P)rini-Q)S{-Qhp, (118) 
7p5(0)r^,„(0)x,^,(P)7., (119) 



and 



L'q^AP) = IpSiQ) K{P) S{-Qhp ,i = l,2. (120) 

NB. It is apparent from the first lines in each of 
Eqs. (116), (117) that there is no mixing between colour 
antitriplet and colour sextet diquarks. 

In this sequence of equations one has arrived at an ana- 
logue of Eq. (80). Although somewhat more complicated 
because of the colour-antitriplet nature of the diquarks, 
the series expresses a recursion relation that may be writ- 
ten compactly as 



An 
A" 



[Tg 



i^l,n— 1 
fi2,n-l 



+ [Tl 



Tl,n-1 
^qqv 

7-2, n-1 
^qqv 



with 



[Tg 



"1-3' 




"3 r 




, [Tl] = - 


1 3 


3 -1 



(121) 



(122) 



and which is subject to the initial conditions: A° = 0. 
NB. In this case, owing to modifications associated with 
the colour algebra, Ap+^(P) 7^ 0. Hence both the scalar 
and axial-vector diquark Bethe-Salpeter equations con- 
tain nonplanar contributions generated by the analogue 
of Eq. (102). 

From this point one proceeds for both the scalar and 
axial-vector diquarks as one did from Eq. (98) for the p- 
meson. The matrices Ai^ and A2p are expressed in terms 
of independent vectors (i^, (32] viz.. 



where A^^ and Ajj^ are Dirac matrices alone. It follows 



Pa 



M^Ai^ , 



(123) 
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FIG. 12: Evolution of the reciprocal diquark masses, calcu- 
lated in the chiral limit, with the parameter C, a positive value 
of which expresses attraction in the colour-octet projection of 
the quark-aritiquark scattering kernel: solid line - scalar di- 
quark, dashed line - axial-vector diquark. The dotted lines 
are explained in connection with Eq. (126). The mass-scale 
g = 0.69 GeV. 



which can be combined into a single column vector P^^. 
That vector is given by 



1 



1-M,TlL, 



■M^TcGgg, (124) 



wherein Gqqi, is the matrix constructed from Eqs. (118) 
and (119) using the fully rcsummcd quark-gluon vertex; 
i.e., T^^ on the r.h.s. in these equations, and Lqqi, is 

defined implicitly via Eq. (120); namely, the right-hand- 
sides are combined into a single matrix and identified as 
Lqqvl^qq- With (3 qq in Eq. (124) one has in hand a closed 
form for 



[2Au 



A 



2^J 



(125) 



in Eq. (107) and therefrom ngqAiP"^) in Eq. (111). Thus 
has one arrived finally at the 2x2 matrices needed to 
construct analogues of Eq. (103); viz., two characteristic 
equations, one for the scalar diquark and another for the 
axial- vector. A bound diquark exists if either or both of 
these equations possess a real zero. 

We have studied the behaviour of these characteristic 
equations as the parameter C in Eq. (22) is varied and 
our results are summarised in Fig. 12. A value of C = 
implements the rainbow-ladder truncation. In this case, 
as we promised, both scalar and axial- vector diquarks are 
bound, with mo+ = 1.55 GeV and mi+ = 1.70 GeV, in 
agreement with Ref. [31]. 

When C is evolved to negative values, which corre- 
sponds to net repulsion in the colour-octet projection 
of the quark-antiquark scattering kernel, the position of 
the zero in both characteristic equations moves deeper 
into the timelike region; i.e., the diquark masses increase. 



This continues briefly until, at C^^ « —0.023 w —1/43, 
the characteristic equation for the axial-vector diquark 
no longer has a solution. For the scalar diquark this hap- 
pens at Cq+ « —0.027 « —1/37. It is therefore clear that 
very little repulsion in the colour-octet quark-antiquark 
scattering kernel is sufficient to prevent the appearance of 
diquark bound states. NB. If one considers single-gluon 
exchange between the quark and antiquark then a value 
of C = —1/8 is obtained, Eq. (19). This is approximately 
five-times larger than these critical values, an observation 
which further elucidates the results in Refs. [31, 32]. 

The evolution of C to positive values provides alto- 
gether new information and insight. To begin, the bound 
diquarks that existed in rainbow-ladder truncation sur- 
vive and their masses decrease continuously with increas- 
ing attraction. Such behaviour, although not necessarily 
welcome, might have been anticipated on the basis of 
continuity, and a decrease in a bound state's mass with 
increasing attraction between its constituents is not un- 
usual. This smooth development continues, for the axial- 
vector diquark until ft: 0.34 and for the scalar diquark 
until Cq+ « 0.36, at which point it changes dramatically. 
The masses suddenly begin to increase rapidly and their 
behaviour thereafter is described by 

^ = 5.59 V0.388 - C , — = 5.75 \/0.363 - C . 



mo+ 



mi+ 



(126) 

Equations (126) are represented by the dotted lines in 
Fig. 12. 

The value of C at which the behaviour of the masses 
changes qualitatively is correlated with the movement of 
the cusp evident in Figs. 1 7 into the domain that affects 
the position of the zero in the characteristic polynomials. 
This means that, with net attraction in the colour-octet 
quark-antiquark scattering channel, the expulsion of di- 
quarks from the bound state spectrum follows immedi- 
ately upon the active expression of confinement by the 
dressed-quark propagator in the bound state equation. 

There are no bound diquarks in the spectrum obtained 
with the value of C = C in Eq. (68) suggested by the 
lattice data. 



V. EPILOGUE 

Herein we have explored the character of the dressed- 
quark-gluon vertex and its role in the gap and Bethe- 
Salpeter equations. Our results are relevant to the mech- 
anism and realisation of confinement and dynamical chi- 
ral symmetry breaking, and the formation of bound 
states. 

We employed a simple model for the dresscd-gluon in- 
teraction to build an Ansatz for the quark-gluon ver- 
tex whose diagrammatic content is expressly enumerable. 
The model reduces coupled integral equations to alge- 
braic equations and thereby provides a useful intuitive 
tool. We used this framework to argue that data ob- 
tained in lattice simulations of quenched-QCD indicate 
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the existence of net attraction in the colour-octet projec- 
tion of the quark-antiquark scattering kernel. 

We observed that the presence of such attraction can 
materially affect the uniformity of pointwise convergence 
to solutions of the gap and vertex equations. For ex- 
ample, in the timelike region, the vertex obtained by 
summing an odd number of loop corrections is point- 
wise markedly different from that obtained by summing 
an even number of loops. The two subseries of vertices 
so defined follow a different pointwise path to the com- 
pletely resummed vertex. Furthermore, correlated with 
this effect, we found that solutions of a gap equation de- 
fined using an odd-loop vertex and those of a gap equa- 
tion defined with an even-loop vertex follow a different 
pointwise path to the solution of the exact gap equa- 
tion; viz., the gap equation defined by the fully resummed 
vertex. This entails that the solutions of two gap equa- 
tions that are defined via vertex truncations or vertex 
Ansatze which appear similar at spacelike momenta need 
not yield qualitatively equivalent results for the dressed- 
quark propagator. This is especially true in connection 
with the manifestation of confinement, for which the be- 
haviour of Schwinger functions at timelike momenta is of 
paramount importance. 

Our study showed that the dependence of the dressed- 
quark-gluon vertex on the current-quark mass is weak un- 
til that mass becomes commensurate in magnitude with 
the theory's intrinsic mass-scale. For masses of this mag- 
nitude and above, all vertex dressing is suppressed and 
the dressed vertex is well approximated by the bare ver- 
tex. 

It is critical feature of our study that the diagrammatic 
content of the model we proposed for the vertex is ex- 
plicitly enumerable because this enables the systematic 
construction of quark-antiquark and quark-quark scat- 
tering kernels that ensure the preservation of all Ward- 
Takahashi identities associated with strong interaction 
observables. This guaranteed, in particular, that inde- 
pendent of the number of loop corrections incorporated 
in the dressed-quark-gluon vertex, and thereby in the gap 
equation, the pion was automatically realised as a Gold- 
stone mode in the chiral limit. Such a result is impossible 
if one merely guesses a form for the vertex, no matter how 
sound the motivation. As a consequence we could reli- 
ably explore the impact on the meson spectrum of attrac- 
tion in the colour-octet projection of the quark-antiquark 
scattering kernel. In accordance with intuition, the mass 
of a meson decreases with increasing attraction between 
the constituents. 

We found that the fidelity of an approximate solution 
for a meson's mass increases with the number of loops 
retained in building the vertex. However, a given con- 
sistent truncation need not yield a solution. This fact is 
tied to a difference between the convergence paths fol- 
lowed by the odd-loop vertex series and the even-loop 
series. In addition, we observed that with increasing 
current-quark mass the rainbow-ladder truncation pro- 
vides an ever more reliable estimate of the exact vector 



meson mass; i.e., the mass obtained using the completely 
resummed vertex and the completely consistent Bethe- 
Salpeter equation. For pseudoscalar mesons, this is even 
more true because the rainbow-ladder and exact results 
are forced by the Ward-Takahashi identity to agree in the 
chiral limit. This is practically useful because it means 
that the parameters of a model meant to be employed in 
rainbow-ladder truncation may reliably be fixed by fitting 
to the values of pseudoscalar meson quantities calculated 
in the neighbourhood of the chiral limit. Moreover, both 
in rainbow-ladder truncation and with the complete ver- 
tex and kernel, the splitting between pseudoscalar and 
vector meson masses vanishes as the current-quark mass 
increases. In our complete model calculation this split- 
ting is 130 MeV at the c-quark mass and only 40MeV 
at the 6-quark mass, a pattern which suggests that the 
pseudoscalar partner of the T(15') cannot have a mass as 
low as that currently ascribed to the r]b{lS). 

Our framework permitted a treatment of diquark cor- 
relations on precisely the same footing as mesons. In this 
case we exposed particularly interesting consequences of 
attraction in the colour-octet projection of the quark- 
antiquark scattering kernel. Colour-antitriplet diquark 
correlations form bound states in rainbow-ladder trun- 
cation. We demonstrated, however, that introducing a 
very small amount of repulsion into the kernel eliminates 
these states from the strong interaction spectrum. If, on 
the other hand, one introduces a small amount of colour- 
octet attraction, these diquarks persist as bound states. 
However, as one increases the amount of attraction to- 
ward that indicated by lattice simulations a dramatic 
change occurs. The diquark masses suddenly begin to in- 
crease rapidly when the amount of attraction is sufficient 
to cause the confining nonanalyticity of the dressed-quark 
propagator to enter into the domain which affects binding 
in the colour-antitriplet channel. The growth continues 
and, with little further increase in attraction, diquarks 
vanish from the spectrum. (NB. There are no bound 
diquarks in the spectrum obtained with the amount of 
attraction suggested by lattice data.) This occurs be- 
cause the colour algebra in the antitriplet channel acts to 
screen attraction between the quarks in a diquark correla- 
tion so that, with sufficient strength in the gap equation, 
dressed-quark confinement alone becomes the dominant 
feature and diquarks are thereby expelled from the bound 
state spectrum. 

It would be natural and useful to extend this work to 
systems composed of quarks with unequal current-quark 
masses. 
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